Heat Capacity as A Witness of Entanglement 
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CO ' We demonstrate that the presence of entanglement in macroscopic bodies (e.g. solids) in thermo- 

dynamical equilibrium could be revealed by measuring heat-capacity. The idea is that if the system 
CO . was in a separable state, then for certain Hamiltonians heat capacity would not tend asymptotically 

to zero as the temperature approaches absolute zero. Since this would contradict the third law of 
thermodynamics, one concludes that the system must contain entanglement. The separable bounds 
are obtained by minimalization of the heat capacity over separable states and using its universal low- 
temperature behavior. Our results open up a possibility to use standard experimental techniques of 
solid-state physics - namely, heat capacity measurements - to detect entanglement in macroscopic 
samples. 
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I. INTRODUCTION 



g 

Entanglement is not only a fundamental and curious feature of purely quantum nature^, but it is recognized as a 
;~r physical resource useful in tasks such as quantum computation, quantum cryptography or reduction of communication 
complexity. A state p of N subsystems is entangled if it cannot be prepared by local operations and classical 
communication, i.e. it cannot be written as a convex sum over product states: p = Y] j ^jPj O ... ® Pj\ where the 

factorizable state p^p ® ... <g> Pj N ^ of N individual systems occurs with the weight Wj > w j = -0 m the mixture. 

A particularly interesting question is whether or not microscopic phenomena, such as quantum correlations between 
individual constituents of a macroscopic body (for example, individual spins of a magnetic solid) may affect its 
macroscopic properties. The usual expectation is that non-classical effects in a macroscopic body vanish due to 
interaction of its many degrees of freedom with the environment (decoherence). In order to experimentally test this 
claim, one needs to apply some types of entanglement criteria, such as the Bell inequalities 2 , to macroscopic bodies. 
Due to limited access to the state of these systems this is, however, usually not possible. Recently, it was shown that 
some thermodynamical properties, such as the internal energy^ii^i, or the magnetic susceptibility 8,9 can detect 
entanglement between microscopic constituents of the solid. They can hence be used as entanglement witnesses^. 
Their additional advantage is their extendibility, that is their proportionality to the size of the sample. For this 
reason we do not need to know how much of the material is studied, but express the quantities as specific (molar, 
per site, etc.). However, the two macroscopic quantities mentioned above also have some drawbacks as entanglement 
witnesses. The magnetic susceptibility can be applied only to magnetic systems and a specific class of Hamiltonians 
(isotropic 9 ). Determining the internal energy at a given temperature might be a complicated experimental task. 

In this paper we show how entanglement in macroscopic samples and in thermodynamical equilibrium can be 
detected by measuring heat capacity. Our method of entanglement detection is both simple and generic. Unlike 
internal energy, measuring specific heat of a solid is a well-established experimental routine in solid state physics. 
Furthermore, heat capacity is a generic property of materials and can thus also be measured on non-magnetic systems 
(in contrast to magnetic susceptibility). 
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In a more general context, our result shows a new link between two fundamental theories, quantum mechanics and 
thermodynamics (see Refs. 11-13 for other interesting links between the two theories). It is related to the Nernst's 
theorem, also known as the third law of thermodynamics^. In the original version, the theorem states that the entropy 
at the absolute zero temperature is dependent only on the degeneracy of the ground state. Alternatively, it can be 
expressed as a requirement of unattainability of the absolute zero temperature in a finite number of operations^. 
This requires the specific heat to tend asymptotically to zero as the temperature approaches the absolute zero. We 
will show, however, that, for certain Hamiltonians and under the assumption that the system is in a separable state, 
one obtains a non-vanishing value for heat capacity (separable bound) as the temperature approaches the absolute 
zero. Since this contradicts the third law of thermodynamics one concludes that the system must be in an entangled 
state. The separable bound for heat capacity is obtained in two different ways, by direct minimization of the value of 
heat capacity over separable states (the explicit example considered is the Ising model in a transverse magnetic field), 
and by referring to the universal behavior of heat capacity close to absolute zero. Using these methods we obtain 
the range of physical parameters (critical temperature and strength of the magnetic field) for which entanglement is 
present in various classes of systems. 

One might question the relevance of our method for entanglement detection since the method requires knowledge of 
the Hamiltonian of the system, and thus one could directly determine its eigenstates, build thermal states therefrom 
and check their separability. As we know, computation of eigenstates is in general a hard problem and the origin 
of some major difficulties in solid state physics. Furthermore, even if the thermal (mixed) state is known, it is in 
general hard to find out whether it is separable or not. Our method requires only knowledge of eigenvalues (partition 
function) and can easily be experimentally implemented. 

Consider a system descried by a Hamiltonian H to be in a thermodynamical equilibrium at a given temperature 
T, Its thermal state is given by pt = exp (— -M^) /Z, where k is the Boltzmann constant and Z = Tr[exp (— -gf)] is 
the partition function. The knowledge of the partition function allows us to derive all thermodynamical quantities. 
For example, the internal energy is given by U — —dlnZ/dp — Tr[H exp (— Mp\\jZ = Tt[ptH], where (3 — 1/kT. 
Similarly, the heat capacity, C = dU/dT — l/(kT 2 )d 2 \nZ/d(3 2 , is proportional to the variance of the Hamiltonian, 

C = ^ = ^((H 2 ) - (H) 2 ). (1) 

We first consider the particular case of an Ising chain in a transverse magnetic field and then consider a more 
general case. 



II. TRANSVERSE ISING MODEL 



A. Introduction 



The Hamiltonian of an Ising ring of N spins- \ in a transverse magnetic field is given by 

N N 

H Ismg =JJ2 °f °f+i + B of- (2) 

where we assume the periodic boundary conditions N + i = i. Here B is the external transverse magnetic field, and 
J denotes the coupling constant, taken in this section to be 1. The Pauli matrices a\ (j — x, y, z) have the following 
actions on ith qubit: crf\ki) = (— l) ki \ki) and af = \ki 1) (k = 0, 1), and © denotes the summation modulo 2. 

In the following we will prove that no product state belongs to the eigenbasis of Hi S i ng . Therefore, the variance 
of Hi s i n g and consequently heat capacity, cannot vanish within the set of these states. Since taking a convex sum 
over product states can only increase the variance, we will conclude that heat capacity cannot vanish for the set of 
all separable states. 



B. Proof of non-separability of the eigenstates 



The proof that all eigenstates of the Hamlitonian Hi S i ng are entangled is by reductio ad absurdum. We first assume 
the opposite, i.e. that at least one of the eigenstates is a product state, and then arrive at a contradiction. Suppose 
that the product state \tp) — <S>fli(ai|0i) + &i|lj)) is an eigenstate that is associated to the eigenvalue E and <ij ^ 
and bi ^ for all i (If this was not the case, i.e. di = or hi = for some i, would not be an eigenstate, because 
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the magnetic field term in H] S i ng flips the state |0) to |1) and vice versa for every qubit). Under the assumption that 
\ip) was an eigenstate with an eigenvalue E, the following would need to hold: 



E = 



(00...0\H Ismg \ij) _ (10...0\H Ising \il>) 



(oo...o|V) 



(io...o|v) 



(3) 



The proof is as follows. The two denominators are equal to Yii=i a i ano - ^7 Y[i=i a i- Under the action of Hising onto 
in the expression on the left, only the terms for which either all spins are in the state |0) or for which one spin 
is in the |1) state and the rest in |0) remain. The numerator on the left is hence equal to Ilili a i 0^ + ^J2iLi ^ 
Under the action of Hi S i ng onto \ip) in the expression on the right only the following terms give a contribution: |10...0) 
(equal to N — 4), |0...0) and the states in which the first and one other spin are anti- aligned to the rest. The second 



numerator results in YiiLi a i {^(-^ 
reduces to a quadratic equation 



4)+B 



1 + £7 (Ei=i 7-- It ) }• After simplifying both sides of Eq. ©, it 



uih. 

at 



4* 
ai 



B = 0. 



(4) 



Note that Eq. Q has two solutions, some ^ — xq and — Interchanging all bras (0| and (1| in Eq. we arrive 

to the same form of the equation as Eq. ((4]), but for the inverse of the fraction, This means that also — x$ and 
must satisfy the equation. This would be true only if xq = ±1, which is, however, not possible. QED. 



C. Minimization of the variance of Hamiltonian over separable states 

We minimize the variance A 2 (Hj S i ng ) over all separable states. To find the minimal value of the variance for all 
separable states it is sufficient to perform minimization over pure product states only. As shown by Hofmann and 
Takeuchii&, any convex mixture over product states p = J2i Wi\4>i) {4>i\ (wi > 0, J2i w i — 1 are weights of the product 
states in the mixture) can only increase the variance of an observable A, with respect to the variances A 2 (A)i of A 
for individual states \<fri): 



A 2 (A) = 5>*<(^-<^» 2 >i = E'' 



(A 2 ) t -(A) 2 + ((A) t -(A) ) 2 \ > ^2 w l A 2 (A) i , 

\ =A 2 (A), >0 



with (A)i — ((f>i\A\(j)i) . This implies that the bound that is obtained for pure product states will also be the bound 
for all separable (in general, mixed) states. 

It can be shown that if the state is a product state, then one has 

(N N N \ 

- ]>>* 2 > 2 <<fi> 2 +25>f)< < 7f +2 ) - 5>f)« +1 ) 2 (af +2 ) 
1=1 2 = 1 2=1 / 

N / N \ 

-ZB^iafKotuKiof) + (<rf +1 )) + B 2 (n — ^(af ) . 

2 = 1 \ 2=1 / 



(5) 



Due to the translational symmetry of the system we expect that the product state, which minimizes this expression, 
is (quasi-)translation invariant. We first consider product states with a period of two sites. In such a state every odd 
spin is in the same state l^i) and, similarly, every even spin is in the state \ip2}- Therefore, the state of N spins (./V 
is here taken even) is \ip) = (|"0i) IV^))* 8 ^ 2 and herein called two-translation invariant. The assumed two-translation 
invariance of the state allows neighboring spins to have anti-parallel z components - to minimize the interaction 
energy - and at the same time to have the x components, all anti-aligned to the magnetic field. Since in Eq. (J5]) one 
has terms dependent on the correlations between non-neighboring spins, one could expect that the variance should 
be minimized over states with a period of more than two sites. We have considered four-translation invariant states, 
|^) = (|V'i)IV'2)|V'3)l'04)) (8A ^ 4 (thus N is divisible by 4) and have shown numerically, that the same bound for the 
variance is obtained as in the case of two-translational invariant states. We assume that the Bloch vectors of j^i) and 
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FIG. 1: The minimal variance per site A 2 (Hi B i ng )/N of the Hamiltonian of a transverse Ising ring versus the magnetic field B 
(J — k — 1). The minimization of the variance is performed over 2-translation invariant product states (l^i)!^))®^ 2 - 
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FIG. 2: (Color online) The specific heat per spin (in units |J| = k = 1) of a transverse antiferromagnetic Ising chain versus 
the temperature T and the magnetic field B^. The values below the red line indicate the region of the (T, B)-space, where 
entanglement exists. There the variance of the Hamiltonian is lower than for any separable state. 



\ip 2 ) lie in the xz-plane, that is Vi {il)i\a y \ipi) = 0. Under these assumptions the variance takes a form of 

A 2 (H Ismg ) = N(l + £ + 4-344) 
- 2BN[z 1 z 2 (x 1 +x 2 )} 

+ B 2 ^(2-x\-xl), (6) 

where we have adopted the notation Xi — (ipi\a x \ipi) = sin^ and Zi = (ipi\a z — cos^. Since the expression |6|) is 
proportional to N, it is convenient to discuss the specific heat capacity, that is the heat capacity per spin. 



D. Results 



Figure I presents the results of a numerical minimization of Eq. ^ over 9\ and 9 2 as a function of the magnetic 
field with units of J = k = 1. The plot confirms that no eigenstate of the Hamiltonian has the proposed form for any 
B 7^ 0. In the limit of a very strong field, all spins tend to orient themselves toward the field and build a product 
state, however, the interaction term contributes the variance with a constant magnitude. This explains why the curve 
in Fig. 1 does not tend to zero as B increases to infinity, but saturates at B w 3.5. 

The results of the optimization are compared to the values for specific heat of an infinite Ising ringii. Katsuraii 
obtained analytical forms of thermodynamical quantities by an exact solution of the Hamiltonian eigcnproblemii. 
The solution was achieved with the Jordan- Wigner transformation followed by the Fourier and the Bogoliubov trans- 



FIG. 3: The specific heat per spin (in units |J| = fc = l)ofa transverse antiferromagnetic Ising chain versus the temperature 
T for B = 2. The gray line, 0.4197/T 2 , is our entanglement witness. 



formations. The heat capacity per spin was found to be 

Ci S i n g 1 f n f(B,w) 



N nT 2 J C osh 2 ^ 



-du, (7) 



where the Boltzmann constant is k = 1 and f(B,ui) = \J\ — 1B cos ui + B 2 . Figure 2 presents the expression ([7|) 
for heat capacity per spin as a function of the external magnetic field and temperature. Values below the red line 
cannot be explained without entanglement; in this region the A 2 (Hj S i ng ) is lower than for any separable state. An 
interesting observation is that with increasing the strength of the magnetic field, the critical temperature below which 
entanglement is detected increases as well. 

Figure 3 presents the heat capacity per spin for a given magnetic field B = 2. The gray line represents a/T 2 , where 
a = 0.4197 is the minimal Hamiltonian variance over separable states for this value of the field. For the temperatures 
below the value of intersection of the two lines the state is entangled. 



III. UNIVERSAL LOW-TEMPERATURE BEHAVIOR AND ENTANGLEMENT 



In the following we will show that the specific heat is an entanglement witness, whenever the internal energy is. 
The argument is based on two results of statistical mechanics, universality and scalability, which together state that 
the partition functions of various systems are similar and characterized by a small number of universal parameters, 
such as the central charge or the dimensionality of a lattice. 



A. Internal energy as entanglement witness 

For the sake of the further discussion, we recall that for certain Hamiltonians internal energy is an entanglement 
witness. This is nicely illustrated in the case of the Heisenberg antiferromagnetic ring Hamiltonian^i, given by 
H X xx — J Si=i ^ ' ( w itb denoting the z-th spin vector of the magnitude s in a ring, periodicity guaranteed 
by i = i + N, and J > 0). Under optimization over product states \ip) — the lowest possible energy is 

Eb = —JNs 2 , since for product states one has: 

N 

E = (iP\H xxx \ip)=jJ2(?)(? +1 ) 

8=1 

JV 

> -jJ2^i^V(^ +1 } 2 > -JNs 2 = E B . (8) 
i=l 

Here, (?) = <^|s*|^). 

By convexity, the same bound holds for all separable states. On the other hand, in the limit of large N, the ground- 
state energy per spin, for example, for s = \ is equal to -0.443J 18 . Thus, in this case for all temperatures below a 
critical one (Tc), for which the internal energy per spin is U xxx .t c = —0.25 J, the state must contain entanglement. 
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B. Gapless systems 



We demonstrate, how the separable bound on the internal energy can be used to derive the separable bound on the 
heat capacity. First, let us consider a class of materials, in which the lowest part of the energetic spectrum is continuous. 
Infinite half-odd integer spin chains and ring are examples of such systems. For sufficiently low temperatures, their 
internal energy (hereafter, "per spin"') can be expanded to a polynome: U(T) = Eq + (/cT) 7 (co + c\kT + ...), with 
E being the ground-state energy, and 7, Co, Ci, ... are material-dependent constants. At sufficiently low temperatures, 
when the higher order terms are negligible, the specific heat is proportional to a power of T as given by: 

C{T) = QT = jcqWIP = 7 — - — . (9) 

Now we use the fact that if the internal energy is bounded for all separable states, i.e., has a separable bound Eb, 
then the specific heat is bounded for all separable states as well. Namely, one has 

C(T)> 7 ^^. (10) 

This allows us to use the heat capacity as entanglement witness. Given a Hamiltonian one determines the ground-state 
energy Eq and the minimal energy over all separable states, Eb- If the ground state is separable, i.e. Eq = Eb, the 
inequality (JTDJ) is in agreement with the third law of thermodynamics, as one can have C — > with T — ► 0. In the case 
Eq 7^ Eb, one compares j Eb ~ e ° with the real temperature dependence of the heat capacity C'(T), obtained either 
from an experiment or from a theory. Since 7 ~ Eb ~" E ° diverges as the temperature approaches the absolute zero and 
C(T) has to tend to zero to be in an agreement with the third law of thermodynamics, the two curves intersect each 
other. The intersection point then defines the critical temperature below which entanglement exists in the system. 
The method does not require the knowledge of the energy eigenstates, but only of the ground state energy and a 
separable bound for the internal energy. 

We note that another thermodynamical entanglement witness has similar features to heat capacity. Namely, the 
magnetic susceptibity also might diverge in the limit of infinitely low temperatures for all separable states^. In both 
cases quantum entanglement is decisive for finiteness of the low-temperature values of the thermodynamical quantities. 

Consider gapless models in 1+1 dimensions. For these systems, the conformal field theor y 19 i 20 i 21 predicts C = 
(nck 2 T) / (ihv) for the low-temperature behavior of the specific heat (throughout this paragraph, we assume J = k = 
h = 1). Here c stands for the central charge of the corresponding Virasoro algebra and v is the spin velocity. For 
half-integer spins c = 6s/(2 + 2s). An example is the infinite xxx antiferromagnetic spin-i Heisenberg chain. From 
inequality (|10[> one obtains that the values of the specific heat below C < 0.386/T manifest entanglement in the 
Heisenberg chain. The approximation C « (2/3)T is valid for temperatures below 0.1—. 

Another example is the xx spin-^ antiferromagnets whose solution is also given in Ref. 17. The internal energy per 
spin of an infinite xx ferromagnet described by the Hamiltonian H xx = J^2i{o x <J x +l + a y a y +1 ) is given by 

4 j r^ 2 ( 2 J \ 

Uxx = / tanh —— sino; I smcoduj. (11) 



tt n V kT 



We will use the following approximation: at sufficiently low temperatures T << k/J, the argument tanh[(2 Jx) /(kT)] 
of the integral is a very steep function, which is well approximated by f(x) — (fy) x — ( j^) x 2 for < x < ^j- and 
J 



f(x) = 1 for x > Thus the internal energy can be written as 



Uxx ~ / f(Lu)ujduj / smtvduj 



-f 

IT JO 


.1 

f(u>)ujduj 


4J 

TT 


5k 2 T 2 




k 2 T 2 


3ttJ 


2J 2 


4J 
h 


k 2 T 2 




TT 


3ttJ ' 





kT 
J 



(12) 

where we have used expansions sin a; w x and cosx m 1 — x 2 /2. From Eq. (|12|) it is clear that the ground-state 
energy per spin is —4 J /it, while the energy bound for separable states is — J. Hence, by Eq. (p~0|) the thermal state is 
entangled if it satisfies C < 2J(4/tt— 1)/T as long as the specific heat can be approximated by a linear function of T. 
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C. Gapped systems 



The integer spin chains have an energy gap A between the ground state and the first excited state, even in the 
thermodynamical limit. The gap is also a feature of systems described by a Hilbert space of finite dimensions. For 
all gapped systems, the low-temperature behavior of the specific heat is given by 

C = c'T s e-zr, (13) 

where d and 8 are some material-dependent constants^ 3 -. The internal energy can then be written as U(T) = 
U a + J T C{T')dT' = U + d(A/k) s+1 T [-6 - 1, A/(fcT)]. Since for large x we have T(a, x) « e~ x x a - 1 ^, we obtain 
U(T) E + dkT s+2 e- A / {kT y& at sufficiently low temperatures. Finally, using the energy bound Eb for separable 
states we obtain the corresponding bound for the specific heat: 

C>^§^ (14, 

for all separable states. Again, one can invoke the third law of thermodynamics to argue for the necessity of existence 
of entanglement at sufficiently low temperatures whenever Eg ^ E . An exemplary gapped system is an infinite xxx 
spin-1 Heisenberg chain with the values d = A 5 / 2 /\/2^, 6 = -3/2, A = 0.411J, and E = - 1.401 j2Si2£, By Eq. ©, 
Eb = — J- Thus, within the approximation range, all low-temperature values of the specific heat below 0.165 J / (kT 2 ) 
cannot be explained without entanglement. 



IV. CONCLUSIONS 



In summary, we have shown that the low-temperature behavior of the specific heat can reveal the presence of 
entanglement in bulk bodies in the thermodynamical equilibrium. For certain Hamiltonians and under the assumption 
of having separable states only, the specific heat would diverge at temperature approaching the absolute zero. This 
might be because none of the eigcnstates of the Hamiltonian are a product state and hence its variance cannot vanish 
within the set of these states or, more generally, because at least the ground state is entangled. In the latter case 
we involve the separability bound and the universal low-temperature behavior of internal energy to argue for non- 
classicality of a thermal mixture in certain systems. One may therefore say that in these systems the validity of the 
third law of thermodynamics relies on quantum entanglement. 

Thermal entanglement of bulk solids might play an important role in the emerging quantum information technology, 
where non-classical correlations were recognized as one of its main resources^. Our method enables to detect 
entanglement using one of the standard techniques in solid-state physics - measurements of the heat capacity. 
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